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g : ABSTRACT 
O 

^ ' Magnetospheres of neutron stars are anchored in the rigid crust and can be 

twisted by sudden crustal motions ("star quakes"). The twisted magnetosphere 
does not remain static and gradually untwists, dissipating magnetic energy and 
producing radiation. The equation describing this evolution is derived, and its 
^ ! solutions are presented. Two distinct regions coexist in untwisting magneto- 

spheres: a potential region where V x B = ("cavity") and a current- carrying 
bundle of field lines with V x B 7^ ("j-bundle"). The cavity has a sharp 



in 



o 



c/5 \ boundary, which expands with time and eventually erases all of the twist. In this 

process, the electric current of the j-bundle is sucked into the star. Observational 

^ \ appearance of the untwisting process is discussed. A hot spot forms at the foot- 

^ \ prints of the j-bundle. The spot shrinks with time toward the magnetic dipole 

\ axis, and its luminosity and temperature gradually decrease. As the j-bundle 

^ \ shrinks, the amplitude of its twist ip can grow to the maximum possible value 

^ \ "ipraax ~ 1- The stroug twist near the dipole axis increases the spindown rate of 

the star and can generate a broad beam of radio emission. The model explains 
OO 

Q I the puzzling behavior of magnetar XTE J1810-197 — a canonical example of 

magnetospheric evolution following a starquake. We also discuss implications for 
other magnetars. The untwisting theory suggests that the nonthermal radiation 



>< 



^ I of magnetars is preferentially generated on a bundle of extended closed field lines 

near the dipole axis. 

Subject headings: plasmas — stars: magnetic fields, neutron 

1. Introduction 

Neutron stars are highly conducting and strongly magnetized. Their extended magneto- 
spheres are anchored deep in the rigid crust and corotate with the star. The magnetosphere is 
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usually assumed to be static in the co-rotating frame or evolving very slowly as the star ages. 
Electric currents are confined to a narrow bundle of open field lines that connect the star to 
its light cylinder (Goldreich & Julian 1969). The main, closed, part of the magnetosphere is 
usually assumed to be current-free and potential, V x B = 0. 

The standard picture of a static potential magnetosphere is reasonable for ordinary 
pulsars, yet apparently it does not describe all neutron stars. In particular, neutron stars 
with ultrastrong fields B ^ 10^^ G (magnetars) are inferred to have dynamic magnetospheres 
(see e.g. reviews by Woods & Thompson 2006; Kaspi 2007; Mereghetti 2008). Their activity 
is believed to be caused by crustal motions relieving internal stressesj^] In contrast to ordinary 
pulsars, these objects show huge temporal variations in luminosity, spectrum, and spindown 
rate. Theoretically, stresses are expected to build up in the deep crust as the star ages 
(e.g. Ruderman 1991; Thompson & Duncan 1995). In particular, the large Ampere forces 
j X B/c inside magnetars can break the crust and shear it in a catastrophic wayl§ Such a 
starquake twists the magnetic field anchored in the crust, creating V x B 7^ and inducing 
electric currents in the closed magnetosphere (Thompson et al. 2000). The currents are 
approximately force-free and fiow along the magnetic field lines, j x B = 0. They emerge 
from the deep crust sheared in the starquake. 

Twisted force-free magnetospheric configurations were studied extensively in the context 
of the solar corona, and these models can be applied to neutron stars. A simple example is the 
self-similarly twisted dipole. It was constructed by Wolfson (1995) and applied to neutron 
stars by Thompson, Lyutikov, & Kulkarni (2002). This and similar force-free configura- 
tions are magneto static solutions. A sequence of such configurations may be constructed by 
changing their boundary conditions, i.e. displacing the footpoints of the magnetic field lines. 
If the footpoints freeze, the configuration freezes as well. At a first glance, this seems to 
suggest that the implanted twist must freeze when the starquake ends, and wait for another 
starquake. 

In fact, the magnetosphere must evolve after the starquake, even though it remains 
anchored in the motionless deep crust. Indeed, energy is continually dissipated in the twisted 
magnetosphere, because the twist current j = (c/47r)V x B is maintained by a voltage 
$e 7^ established along the magnetic field lines. Thompson et al. (2000) estimated voltage 
$e assuming that the currents are carried by electrons and ions lifted from the star's surface 
against gravity. Beloborodov & Thompson (2007; hereafter BT07) found that $e is regulated 



^Well-studied ordinary pulsars (Crab, Vela) also show glitches in their spindown rates, which are associ- 
ated with sudden crustal deformations. 

^In some cases, the crust may move plastically. 
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by an dischargeo The voltage is significant — comparable to 1 GeV — and implies a 
modest lifetime of the twist, comparable to one year. 

The untwisting dynamics of the magnetosphere remained, however, unknown. Usually, 
resistivity in a plasma leads to diffusion of currents across the magnetic field. Voltage $e 7^ 
imphes an effective resistivity, and one could expect the decaying twist to spread diffusively 
across the magnetosphere. This expectation is incorrect as will be shown below. 

The goal of this paper is to develop an electrodynamic theory of twisted magnetospheres 
that describes their evolution. We focus on axially symmetric configurations. In this case, 
the twist is created through a latitude-dependent azimuthal rotation of the crust. An intro- 
ductory description of twisted magnetic configurations is given in § 2, and their untwisting 
dynamics is qualitatively discussed in § 3. In § 4 we derive the electrodynamic equation for 
axisymmetric magnetospheres. 

§ 5 presents solutions to the evolution equation and explores the mechanism of untwist- 
ing. Observational effects of this process are described in § 6. § 7 compares the theory with 
the recent observations of a starquake in the anomalous X-ray pulsar XTE J1810-197. § 8 
summarizes the results of the paper and discusses implications for magnetars. 



In spherical coordinates r, ^, 0, magnetic field can be written as the sum of poloidal and 
toroidal components. 



where e^, e^, are unit vectors pointing in the r, 6*, directions. We assume that the 
magnetic field is symmetric about the polar axis, i.e. B does not depend on 0. The 
axisymmetric field can be viewed as a foliation of magnetic flux surfaces. (Each surface 
may be obtained by rotating a field line around the axis of symmetry.) 

Let / be the magnetic flux through a circular contour of fixed r = const and 6 = const. 
The function f{r,6) is constant on a flux surface. As is usual in plasma physics, we will 
use / to label the flux surfaces. Note that flux surfaces extending farther from the star 
have smaller /, and / = corresponds to the polar axis 6 = 0. We focus on the closed 



2. Twisted axisymmetric magnetosphere 



B = Bp + = BrBr + Boeo + B^e^, 




The discharge on closed field lines differs from that on open field lines (see Arons 2008 and Beloborodov 
2008 for a recent discussion of of the polar-cap discharge in ordinary pulsars). In both cases, however, the 
discharge is ultimately driven by the magnetic twist V x B 7^ that imposes an electric current. 
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magnetosphere in this paper and neglect the narrow bundles of open magnetic field lines; 
effectively, rotation of the star is neglected. 

Let R be the radius of the star. The magnetic field is force-free outside the starQ 
j X B ^ aX r > R. If we follow the magnetic field lines into the star, significant deviations 
from the force-free condition appear. The 0-component of j x B remains, however, small 
as the crust is relatively fragile to axisymmetric azimuthal displacements (which involve no 
compression). This component can be written as x Bp where jp is the poloidal component 
of the current density. Let Tc be the radius of the lower crust that is strong enough to sustain 
significant azimuthal Ampere forces jp x Bp/c 7^ 0. Outside this radius we assume 

jp X Bp ^ 0, r > r^. (2) 

A rough estimate Vc ~ 0.9R is sufficient for the purposes of this paper (the exact Vc depends 
on the strength of the magnetic field). The effective footpoints of the magnetospheric field 
lines sit at r = r^. We focus in this paper on the region r > and call it "force-free" (in 
the restricted sense jp x Bp = 0). 

The spheres r = R and r = rc define two special fiux surfaces (Fig. 1): 

1. Flux surface touches the surface of the star|§ Flux surfaces f > fn are confined to 
the star, and flux surfaces / < extend outside the star and form the magnetosphere. fn 
represents the total magnetic flux emerging from the star in the region of positive polarity 
(where Br > 0). 

2. Flux surface fc touches the boundary of the inner crust r = r^. Flux surfaces f > fc are 
conflned to the inner crust and the core of the star. 

A twist is pumped into the force-free region r > rc when the footpoints of fleld lines at 
r = Vc are displaced by a starquake. The crust is practically incompressible, and hence any 
axisymmetric starquake produces a pure azimuthal displacement (an axisymmetric displace- 
ment in the 6'-direction would imply compression). 



^Electric fields maintaining the currents in a twisted magnetosphere are relatively weak and do not spoil 
this approximation. 

^ Fig. 1 assumes for simplicity that there is only one such flux surface. This is the case for a bipolar 
magnetosphere, with two regions on the star's surface with opposite polarities of the magnetic field (i.e. 
opposite signs of B,.). One can imagine more complicated axisymmetric magnetic configurations with many 
rings of opposite polarities on the star's surface. The development of electrodynamic theory would be similar 
in those cases. When calculating the evolution of currents on a given flux surface / we would need to know 
the magnetic field only between the two nearest flux surfaces that touch the surface of the star. 
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Fig. 1. — Poloidal cross section of an axisymmetric magnetic configuration. The star is the 
sphere of radius R {shaded). The lower crust is inside radius {dashed circle). The figure 
shows the magnetic axis / = 0, flux surface that touches the surface of the star, and flux 
surface /c that touches the sphere Tc. The magnetosphere is composed of nested closed flux 
surfaces f < Jr] one such flux surface is shown in the figure. 
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Let us call the footpoints of positive polarity [B^ > 0) northern. Consider an initially 
pure poloidal magnetosphere and suppose a starquake shifts the northern footpoints of mag- 
netic field lines through angle A0„(/) and southern footpoints through A0s(/). The created 
twist is described by the relative angular displacement ip{f) = A(/)s — A0„ = — This 
angle can be expressed as an integral along the closed field line: = 0s — 0„ is accumu- 
lated as we move along the field line from its northern footpoint at r = rc to the southern 
footpoint. An infinitesimal displacement dl along the field line corresponds to azimuthal 
displacement h^dcp = (B^/B) dl where = rsin^. Therefore, the twist angle is given by 



I ^ lr>rc B h^' ^ 



where the integral is taken along the field line outside r^. Creation of z/' 7^ implies the 
appearance of toroidal magnetic field B^ in the magnetosphere. 

The toroidal field Brf){r,6) determines the circulation of B along the circular contour 
r = const, 6 = const. By Stokes' theorem, it is related to the electric current / flowing 
through the contour]^ 

= (4) 
cr sm t7 

/ is determined by the poloidal component of the current j^, and / is determined by the 
poloidal component of the magnetic field Bp. At r > the condition x Bp = implies 
that the poloidal currents flow along the poloidal flux surfaces. Therefore, J is a function 
of /. Note also that the definition of I is similar to that of / except that B is replaced j, 
which implies 

dl j 

Any axisymmetric force-free field outside the star satisfies the Grad-Shafranov equation 
that expresses the condition B x (V x B) = in terms of / and /. Its exact solutions 
(matching an interior non-force-free solution) are needed to describe twists with large angles 
ip. Note that configurations with ^ 1 are not expected as they are unstable (Uzdensky 
2002 and refs. therein). If the twist grows beyond the instability threshold ipma^ = C>{1), 
the magnetosphere becomes kink-unstable and ejects a closed plasmoid, which prevents the 
twist growth above ^/'max- 



^The current is maintained through a magnetospheric discharge which fluctuates on a very short (hght- 
crossing) timescale r/c (BT07). We consider the time-average / and treat it as a quasi-steady current. Its 
evolution caused by resistivity is slow (year timescale) and can be viewed as a slow progression through a 
sequence of steady states. The displacement current vanishes in a steady state, so V x B = (47r/c) j. 
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For moderate twists -0 < 1 one does not have to solve the Grad-Shafranov equation. 
Instead, a simple linear approximation may be sufficient: the configuration can be thought 
of as a linear superposition of an initial non-twisted poloidal field Bq and a toroidal field 
< Bq that was created by the footpoint displacement ip. The appearance of B^f, does not 
affect the poloidal field in the linear order — the poloidal correction to B is quadratic in 
Brf,/B. The linear approximation may break at large distances from the star (Low 1986), 
however it describes well most of the magnetosphere. Wolfson & Low (1992) found that the 
relation between I and / obtained in the linear approximation works well even for twists 

■0 ~ 1. 

During the starquake, the twisting motion of the footpoints at r = Tc pumps energy 
into the magnetosphere, which may be released later. The free energy of twisted force-free 
configurations was extensively studied (e.g. Aly 1984). For linear twists, B^ < B, the free 
energy simply equals the energy of the toroidal field component B^, 



The current / through the magnetosphere is maintained by electric field E\\ ^ (parallel 
to B), which implies Ohmic dissipation of the twist energy, E- j 7^ 0. Thus, the magnetic field 
must be gradually untwisted (even though it remains anchored in the static deep crust), and 
eventually the magnetospheric current must vanish. On the other hand, the ejected current 
cannot disappear because it emerges from a static and almost ideal conductor — the deep 
crust r < Tc, where the magnetic field and electric currents remain unchanged. We conclude 
that Ohmic dissipation must re-direct the ejected poloidal current so that it closes below the 
surface of the star. The current must be re-directed across the magnetic flux surfaces, which 
can happen only in the transition layer r Tc between the heavy static conductor and the 
force-free region. Then the current does not penetrate the force-free region on flux surfaces 
f < fn s-nd avoids the magnetospheric dissipation. 

It is instructive to consider an idealized problem where the entire star is a perfect 
conductor, so that E\\ 7^ only outside the star. When the magnetospheric dissipation 
is completed, = on flux surfaces f < Jr in the region r > Vc- The initially ejected 
current now flows in a current sheet inside the star (Fig. 2). The current sheet serves as 
a screen between the twisted field inside the star and the untwisted (potential) field in the 
magnetosphere. When the finite resistivity of the crust is taken into account, the state shown 
in Figure 2 is not final. The current sheet in the non-ideal conductor will acquire a non-zero 




(6) 



3. Fate of the ejected current 
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thickness and gradually spread to the inner flux surfaces until the currents reach deeper crust 
with so high conductivity that it can be treated as an ideal conductor on timescales equal to 
the age of the star. The currents that initially emerged during the starquakc will eventually 
close deep under the surface of the star. The tendency of currents to diffuse toward regions of 
higher conductivity was observed in numerical simulations of Ohmic dissipation in neutron 
star crusts (Sang & Chanmugam 1987). 

In summary, two stages are expected in the evolution of a magnetic twist created by 
a starquake: (1) The current ejected by the starquakc into the magnetosphere is gradually 
drawn into the star. Most of the twist energy is released at this stage. (2) The current 
spreads into deeper layers of the star and eventually collects near the highly conducting 
inner crust. The second, subsurface untwisting stage is much slower because the resistivity 
inside the star is smaller than the effective resistivity of the magnetosphere. We shall focus 
below on the faster flrst stage and treat the star as an ideal conductor. 



4. Twist evolution equation 

4.1. Resistive evolution 

The evolution of magnetic held is related to electric fleld E according to the induction 
equation. 

We will express V x E in curvilinear coordinates deflned as follows. Let us label magnetic 
fleld lines on a flux surface / by the azimuthal angle 0o of their northern footpoints at r = Tc 
(northern footpoints have > 0). Thus, the set of all fleld hues is parameterized by two 
coordinates / and 0o- Let s be a parameter running along the field line (increasing in the 
direction of B). The parameter s may be chosen, e.g. equal to length / measured along the 
field line from its northern footpoint (then \es\ = 1). The coordinate system = (s, /, 0o) 
covers the entire magnetic field that passes through the sphere r = Vc and emerges in the 
force-free region. 



The electric field can be written in components in the new coordinate basis, 

e. = |l. (8) 



i=i ^ 
The length of the basis vectors will be denoted by hi — \ei\. Note that 

dr 



dr 
d(f)o 



f,s 



d(t) 



= e^ = h^e^. (9) 



Fig. 2. — Pattern of electric currents after a global starquake that has twisted the entire 
magnetosphere. Poloidal electric currents are shown by arrows; they flow along the magnetic 
flux surfaces outside Tc (dashed circle). Left: initial twisted state. Right final untwisted 
state. When the magnetospheric dissipation is completed, all currents are sucked into the 
star and close below its surface. Thick curve shows the poloidal cross section of the current 
sheet formed inside the star. It extends from the axis along the sphere r ^ rc, across the flux 
surfaces / < /r. The current sheet turns where it reaches the flux surface /r, and continues 
along this fiux surface. 
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We will need the determinant g of metric gik = ■ e^, which is given by ^ = ■ (ej x e<^). 
Consider an infinitesimal axisymmetric ring perpendicular to the poloidal component of B. 
Its surface element is ej x Sfpdf d(f), and the area of the ring is 27r |e/ x e^\df. The magnetic 
flux through the ring, df, is 

df = 27TB-{efXe^) df, (10) 
which implies the identity 27rB ■ (e^ x e^) = 1. We substitute B = Beg/hs and flnd 

The general expression for V x E in curvilinear coordinates is 

{VxEy=y—^, (12) 

where e^^'' is Levi-Civita symbol and -Ej = ■ E are the covariant components of the electric 
fleld in coordinate system q\ Using = /ijej, we obtain the azimuthal component of V x E 
in the normalized basis e, and flnd 



c dt hg \ df ds 

Let us divide both sides of this equation by 2TTBh^/hs and integrate it over s along an entire 
closed fleld line (including its part at r < Tc). Then the second term on the right-hand side 
disappears, and we get 

1 fdB.h.ds d f^^^^_ ^^^^ 



2nc J dt h^B df 
The integral on the right-hand side is the net voltage induced along the magnetic fleld linejil 

$e = j Egds = j E\\ dl, (15) 

where = E ■ and dl = hgds is the length element along the fleld line. Equation 0141) 
shows that the twist evolution is controlled by the longitudinal voltage, as expected. Using 
B^ = 21{f,t)/ch^ (eq. H]), we rewrite this equation as 

J_ I\dlB^^ . . 

27X0 J I dt Bhs df ' ^ ' 



This voltage is not electrostatic and does not vanish for a closed contour. It is the self-induction voltage 
of the gradually decaying twist, see § 2 in BT07. 
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Note that dB^/dt = and dl /dt = a.t r < Vc since the magnetic field remains static inside 
the static ideally conducting inner region, dl/dt jumps from to its value in the force- free 
region at r ^ r^, in a transition layer of thickness A <C r^. The contribution from this 
layer to the integral on the left-hand side of equation flTB]) is small and we neglect it. Then, 
effectively, the integral is taken only along the force-free part of the field line at r > rc. The 
current / and its time derivative dl/dt are constant along this part of the field line. Then, 
using the expression for the twist angle (eq. [3]) , we find 

27icl dt df ^ ' 

Equation (fT7|) describes the evolution of I{f,t) and the corresponding ipifyt)- The twist is 
changing with time as the magnetic field lines gradually slip in the resistive magnetosphere 
and connect new footpoints with different 0s — 0„. Thus, ip is changing despite the fact that 
the magnetosphere remains anchored in the static deep crust. 



4.2. Linear twists 



For small twists < 1, the magnetic field can be written as B ^ Bo+ B^(t) (§ 2), where 
poloidal field Bp = Bq remains static and < Bq. Then B = Bo[l + 0{B'^/B'^)] ^ const 
and the twist evolution equation (ITBi) simplifies to 

1 dip 9$e , , 



2nc dt df 

This approximate equation quickly becomes accurate for B^j^ ^ B: its error is decreasing as 
{B^/By. In the following sections we will use equation (ITSl) to study the evolution of the 
twist amphtude ip. The exact evolution equation ( fTTl) would have to be used when the effect 
of the twist on Bp is of interest (e.g. for calculations of the spindown rate of the star). 

The linearized description of twisted configurations is useful for the first calculations 
of the resistive untwisting. However, the limitations of this approximation should be kept 
in mind. The linearized description may fail at large distances from the star (Low 1986). 
Besides, the force-free configuration may be unable to smoothly adjust to the growing twist. 
Sudden relaxation to a new topological configuration is possible, with partial opening of the 
field lines and the loss of connectivity between the footpoints. The twist evolution equation 
derived above (linear or nonlinear) does not describe such transitions. 
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4.3. Energy conservation law 

The energy of a linear twist < Bq is given by 

Bl 1 f f f 21 1 



Differentiating with respect to time and using equations (IT7|) and (ITS]) , we get 

dE,. (20) 



dt Jq df 

Integrating by parts and taking into account that /(O) = and $e(/c) = [E^^ = on flux 
surfaces confined to the perfect conductor), we find 



dE, 



tw 



dt 



<^edl. (21) 



The right-hand side represents the net Ohmic losses. i?tw is the free energy of the twist that 
is gradually dissipated as the magnetic field evolves according to equation ([T8|l . 



4.4. Magnetosphere with moving footpoints 

The evolution equation derived above assumes that the magnetosphere is anchored in 
the deep static crust following a starquake. It is straightforward to generalize equation f|T8l) 
for magnetospheres with moving footpoints, 

2-1^ (22) 

where uj = dcps/dt — d(f)n/dt is the differential angular velocity of the northern and southern 
footpoints of the magnetic field lines. A crust that remains motionless at all times except a 
sudden starquake at t = is described by 

coif,t)=Mf)m, (23) 

where 6{t) is the Dirac function and ipo is the amplitude of the twist imparted by the 
starquake. The impulsive twisting is a good approximation for recurring starquakes if the 
time between subsequent starquakes is longer than the timescale of untwisting. In the 
opposite limit, the crust is frequently deformed my mini-starquakes or moves plastically. 
Then the footpoint motion can be described by a continuous function uj{f,t). In this paper, 
we focus on the case of impulsive twisting described by equation fl23|) . 
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4.5. Twisted dipole 

For the study of the untwisting mechanism in the next section, it is useful to consider a 
concrete simple magnetic configuration. We will consider a dipole field with the symmetry 
axis passing through the center of the star. Let jJ be the dipole moment. The poloidal flux 
function for the dipole is given by (Appendix A), 

sin^6' 27111 , 
/(r,^) = 27r/x = -^, (24) 

T -'I'max 

where -Rmax = r/ sin^ 9 is the maximum radius reached by the flux surface passing through 
given r, 6. Hereafter, instead of / or -Rmax we will label flux surfaces by the dimensionless 
coordinate 

«-f = ^. (25) 

The last flux surface in the force- free region (marginally emerging from the inner crust: 
-Rmax = rc) has u = Uc = R/tc ~ 1.1. The region < m < 1 corresponds to the magneto- 
spheric flux surfaces. In this region, u = sin^6'i, where 9i is the polar angle of the northern 
footprint of the fleld line on the star surface r = R. 

Suppose now that the fleld has been twisted by a starquake that was symmetric about 
the dipole axis and resulted in a differential rotation of the crust through angle ip{u) (§ 2). 
The poloidal current I and twist angle ip will be viewed below as functions of u and time t. 
The following relation holds between ip and / 



, . ,26) 



(see Appendix A). The twist evolution equation f|T8l) becomes, 

cR d% 



dt n du 

or 



e 



(27) 



dt m^Jl-u/u, du 
The current density in the twisted dipole magnetosphere is given by (cf. eq. [5]), 



5. Mechanism of untwisting 



This section will explore the mechanism of untwisting for twists created by axisymmetric 
starquakes in the dipole magnetosphere (§ l4.5p . Their evolution is described by equation fl27|) . 
Before this equation can be solved, the voltage $e must be specified. 



Consider first what would happen if the magnetosphere was filled by a medium with a 
fixed conductivity a; we will assume in this toy model a{r) = const. Then the current density 
is related to £"11 by Ohm's law j = aE\\ . The star will be modeled as a perfect conductor, so 
E\\ = is assumed inside the star. Then the voltage $e along a magnetospheric field line is 
given by 



(The integral has been calculated using eqs. [29] and [671) Substitution of this result to the 
twist evolution equation ([28ll yields the equation for I{u,t), 



Given an initial twist with current function I{u, 0), one can calculate the evolution of I{u, t) 
by solving this differential equation with two boundary conditions: /(O) = and /(I) = 
const. The latter condition is valid as long as the perfect-conductor approximation is used 
for the star. 

Equation (15T]) is of diffusion type. It has a special feature: the effective diffusion coeffi- 
cient is proportional to ^/l — u and vanishes at u = 1. This fact allows one to qualitatively 
understand the evolution of I{u,t) before solving the equation numerically. It is instructive 
to consider the analogous problem of particle diffusion with a position-dependent diffusion 
coefficient D{u) that vanishes at the boundary u = 1. The vanishing of D means that 
particles "stick" to the boundary. With time, more and more particles get stuck, and even- 
tually the particle density vanishes everywhere except at the boundary. The magnetospheric 
current behaves in a similar way. Far from the boundary, it simply spreads diffusively: 
dl/dt = const d"^ I / du^ aX u <^ 1. At the same time, near the boundary m = 1, the 
current is sucked toward u = 1. The current keeps accumulating at u = 1 until / = at all 
u <l. 

For illustration, we solved numerically equation f[3T]) for a twist that initially has a 
uniform amplitude V'o('") = 0.2. The corresponding initial current function /(u, 0) and 



5.1. Twist evolution in a medium with fixed conductivity 




(30) 




(31) 
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dl /du oc j/B are shown in Figure 3. The evolution of this twist is shown in Figure 4. The 
characteristic diffusion timescale is ta^ = R'^c^/ina, and we express time in units of t^^. As 
expected, the current tends to spread to smaller u and, at the same time, it is quickly drawn 
into the current sheet at m = 1. The sum of the currents flowing in the magnctopshcrc and 
in the current sheet, /tot, remains constant. Our numerical model assumes V = at u > 1 
(inside the star), which allows the current sheet to persist at u = 1. A real star has a finite 
conductivity and the current sheet will slowly spread into the star (toward larger u > 1). 

Note that the twist amplitude ip grows near u = 0, because d^e/du > in this region. 
A strongly twisted bundle develops near the magnetic axis, however, its thickness shrinks 
with time, so its net current decreases. At late times, the twist growth near the axis enters 
a self-similar regime. The amplitude peak ippea.k would grow indefinitely in the limit t — > oo, 
if the twist remained stable. In fact, the growth must be stopped by the MHD instabihty 
expected at ^ 1. An upper limit on ip is set by the field- line opening, which may give 
a partially opened configuration with a lower energy (Wolfson & Low 1992). Besides, the 
strongly twisted configuration is prone to kink instability, which is likely to reconnect part 
of the twisted region away from the star. These dynamic processes are not described by our 
model. We shall assume that they keep ip below V'max = C'(l)- 



5.2. Threshold voltage 

The toy model discussed in § 5.1 is deficient: the real magnetospheres of neutron stars 
are not filled by a medium of a fixed conductivity. Instead, the magnetospheric currents 
are maintained through a discharge with a threshold voltage that is approximately the same 
for any current density j ^ (or, more precisely, for any j exceeding a small value j^). 
The threshold nature of the magnetospheric voltage changes the evolution of the twist. 
Remarkably, it simplifies the solution of the twist evolution equation: the partial differential 
equation will be reduced to an ordinary differential equation. 

We denote the threshold voltage by V. Once <I>e reaches V, copious particle supply is 
available to carry any large current. The particle supply and voltage regulation in magnetars 
was studied in BT07. In principle, there are two sources of particles: (1) the surface of the 
star and (2) creation in the magnetosphere. The current is carried by charges of both 
signs, since the net charge density must be nearly zero to avoid huge voltages. If no are 
created, electrons and ions must be lifted from the star's surface (if ions are available in an 
atmospheric layer atop the solid crust). Maintaining a fiow of ions along a magnetic loop 
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Fig. 3. — A uniform twist created by a global starquake involving the entire magnetosphere 
< u < 1. The twist amplitude in this example is ipoiu) = const = 0.2. The figure 
shows the corresponding current function Io{u) = I{u, 0) and its derivative dio/du, which is 
proportional to current density j (see eq. [29]). Current / is measured in units of / = cfi/AR"^, 
where fi and R are the magnetic dipole moment and the radius of the star. The coordinate 
u = f / fn = R/Rmax labels fiux surfaces (see eq. [25] and Fig. 1); -Rmax is the maximum 
radius reached by the fiux surface, -u = on the magnetic axis (-Rmax — oo on the axis). 
For magnetospheric fiux surfaces m < 1 and u = sin^ 9i where 9i is the polar angle of the 
northern footprints of the fiux surface on the star. Flux surfaces with u > 1 close inside the 
star; this region is shaded in the figure. 
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Fig. 4. — Evolution of the twist (with the initial configuration shown in Fig. 3), calculated 
under the assumption that the magnetosphere is filled with a medium of a fixed conductivity 
cr, and the star is an ideal conductor. Three panels show the twist at times t/to- = 0.01, 
0.1, and 0.7, where to- = -R^c^/47r(T. Solid curve shows the current distribution dl{u,t)/du; 
I in the figure is measured in units of / = cfi/AR"^. The initial dIo{u)/du is shown by 
dotted curve (from Fig. 3). Shaded area represents the magnetospheric current and the 
current sheet Ji at m = 1. The total current /tot = h + h is conserved. Ji grows with time, 
and the magnetospheric current decreases and spreads toward the magnetic axis u = 0. 
Dashed-dotted curve shows the twist amplitude tp. It grows near the axis. 
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requires a minimum voltage, 



GMrrii GMrrii 

eVei = — , (32) 

it rtn 



"-max 



where i?max is the maximum radius reached by the loop and is the ion mass. This 
voltage is ~ 0.2miC^ for loops with -Rmax ^ R- The plasma is lifted into the loop by the 
self-induction electric field (BT07). The exact solution was obtained for this process in 
the one-dimensional circuit model. It shows that voltage keeps growing even after reaching 
Veil and the circuit evolves toward a global double- layer configuration. This result may 
not hold in the complete 3D model that includes the excitation of transverse waves in the 
magnetosphere. Nevertheless, it suggests that lifting plasma from the surface is not the 
ultimate regulator of the voltage. 

A robust mechanism for limiting the voltage is the discharge. If $e exceeds a certain 
threshold V±, the exponential runaway of pair creation occurs (cf. Fig. 5 in BT07), and the 
pairs screen E\\. In a magnetar magnetosphere, the avalanche is triggered when the 
accelerated electrons resonantly scatter X-rays streaming from the star, and the scattered 
photons convert to off the magnetic field. The threshold for the discharge is given by 

eV± ^ 7resmec^ ^ ^ I ( ( GeV, (33) 

where uJx is the typical frequency of target X-rays and 7res ~ {B / BQ){meC^ / tkvx) is the 
electron Lorentz factor at which the electron begins to scatter > 1 X-rays as it travels along 
the magnetic loop. Here Bq — mlc^/eh 4.4 x 10^^ G. Numerical experiments in BT07 
show that the discharge is intermittent on a timescale ~ r/c and proceeds in the regime 
of self- organized criticality. The time-averaged current equals the current imposed by the 
magnetospheric twist, (c/47r)V x B, and the time-averaged voltage $e is close to V±. 

The discharge voltage remains almost independent of the imposed current j unless j is 
reduced below j^. The value of is unknown but small. It may be comparable to cpcj, where 
Pgj is the corotation charge density (Goldreich & Julian 1969) that should be maintained 
in the magnetosphere in the absence of electric currents. For the purposes of the present 
paper, the following description for $e(i) will be sufficient. 



where is much smaller than the characteristic currents induced by the starquake. As will 
be shown below, neither the value of nor the behavior of $e(j) in the transition region 
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j ~ matters. In essence, $e = V(M)0(j), where G is the Heaviside step function. In 
computer simulations, we use a smoothed step function, 

= Viu) W(^), W{x) = arctan(l/A) + arctan [(a: -1)/A] 

^ ' \jj arctan(l/A) + 7r/2 ' ^ ' 

where A ^ 1. Note that the discharge voltage V can be different for different flux surfaces, 
i.e. V in general depends on u. There is a sharp drop in V{u) at m = 1 (voltage is small inside 
the highly conducting star). In computer simulations, we model this drop by introducing 
the factor exp{[e/(l — with e ^ 1. The exact form of this factor plays no role for the 

twist dynamics. 



5.3. Expanding cavity 

Consider again the twist shown in Figure 3, and let us calculate its evolution with the 
new threshold relation between j and $e (eq. IMI) . The numerical solution of equations fl28|) . 
(12^ . and fl35p is shown in Figure 5 for the simplest model that assumes V(m) = const in the 
magnetosphere and V{u) = inside the star. 

Initially, the drop in $e near u = 1 is very steep, i.e. d^e/du is large and negative, 
and hence the current density here is quickly reduced (cf. eq. [27|) . The reduction continues 
until j ~ j-^, which permits $e < V. Then a smooth profile of < $e(w) < V is established 
in a region u-^{t) < u < 1. The threshold nature of the discharge leads to formation of two 
distinct regions in the magnetosphere: 

(1) "Cavity" u^, < u < 1 where < $e < V and j ~ (essentially zero current). The 
current originally injected in this region is sucked into the star and flows in the current sheet 
at M = 1. 

(2) Region u < u^, where $e = V. Here the current remains equal to its initial value at t = 0, 
i.e. the twist remains static^ 

The voltage profile in the cavity has d^e/du < which implies dl/dt < (eq. Hence 
the cavity must grow, as indeed seen in the simulation. The structure of the untwisting 
magnetosphere is shown in Figure 6. 

The boundary of the cavity forms a sharp front, which resembles a shock wave. The 
front starts as a tiny arc emerging from the star at the magnetic equator, and propagates 



^ This is a consequence of V(u) — const. The twist at u < m^, will not remain static if V(m) 7^ const as 
discussed below. 
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Fig. 5. — Evolution of the magnetospheric twist with initial uniform amplitude tpQ = 0.2 
(Fig. 3). The evolution is caused by the discharge voltage $e that is described by eq. flM|) . 
with V(m) = const at m < 1 and V{u) = at m > 1. Dashed curve shows ^e{u); it vanishes 
in the shaded region u > 1 (inside the highly conducting star). Dotted curve shows the 
initial current distribution dlo/du. Solid curve shows the current distribution 81/ du at 
time t = O.Olty {upper panel) and t = O.OAty {lower panel), where ty = ^/cRV; current 
is plotted in units of J = cfi/AR'^. Immediately following the starquake {t = 0), a cavity 
with j ^ ji, ^ forms at m = 1 and grows with time. Its boundary — the current front — 
moves to smaller u i.e. larger i?max = R/u, erasing the magnetospheric currents. Shaded area 
shows the magnetospheric current I^, and the current sheet /i at u = 1. The total current 
-^tot = I* + h is conserved: the erased magnetospheric current flows in the current sheet. 



Fig. 6. — Structure of untwisting dipole magnetosphere (poloidal cross section). The mag- 
netospheric currents are confined to the j-bundle u < u^, (the outer region shaded in yellow). 
Field hues shown in green are potential (V x B = 0) and form the inner cavity with j 0. 
The cavity is bounded by the current front (located at the magnetic flux surface u = m*, 
shown by red curve). The front expands with time, moving to flux surfaces closer to the 
magnetic axis, and the j-bundle shrinks. The erased currents (that initially flowed in the 
cavity) are closed inside the star and flow in the current sheet {thick blue curve). The final 
state shown in Fig. 2 is achieved when the current front u^, reaches the magnetic axis u = 0. 
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toward smaller u (i.e. outward and poleward, to flux surfaces with larger -Rmax)- We denote 
its instantaneous position by u^(t). The profile of the front — the shape of j{u) near — 
is controlled by the behavior of $e at j ~ However, we are not interested in the exact 
profile of the front. It is sufficient to know that it is steep and can be treated as a step 
function. The quantity of interest is the speed of the front propagation du^,/dt. In the limit 
of small (steep front), one can use 

$e = V(M)e(j), (36) 

and derive an explicit expression for du^,/dt (see Appendix B), 

V(M*j \ = 

du^ _ l-u^ a/1-m*/mc . . 

5 l-J ') 



(? du 



d u^V{u) 

+ 



du J\ - u/uc 



where V' = dV /du and Iq{u) = I{u, 0) is the initial current function. This ordinary differen- 
tial equation can be solved for u^,{t). If V = (as in the model in Fig. 5) the front equation 
simplifies to 

du^ _ c^Vu^jl + 1.4mJ 

dt ~ AR{l-u^){dh/du\J ^ ' 

The history of the front propagation is shown in Figure 7. The front starts with extremely 
high speed near m = 1, then decelerates and approaches the axis m = with du^,/dt = 
—cRV /2iJiipQ. It reaches m = (and erases all of the twist) at tend = ^i'o/cRV. 

The current function / of the evolving twist is given by equation fl75p in Appendix B. 
The corresponding twist amplitude is given by 

cR 

ipo{u) -\ V{u)t 0<u<u^ 

Hu.t) = { ^ (39) 
1 hit) u^ <u <1 



cjjiu^ y u, 

where /*(t) = /['U^(t), t] is the net current that fiows through the magnetosphere at time 
t. Equation (1391) together with M^(t) gives a complete analytical description for untwisting 
magnetospheres (in the linear-twist approximation, see § 4.2). 

Note that the numerical simulation shown in Figure 5 gives $e/V < 1 for n ^ 1, i.e. for 
fiux surfaces with footpoints near the magnetic axis. This is not predicted by the analytical 



^ The profile of the front is controlled by the form of function W{j/ji,) in eq. (P5|): in the simulation 
shown in Fig. 5 we chose j^, — lO^^cSo/Svri? and A — 0.2. 
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Fig. 7. — Lower panel: propagation of the current front Ui,{t) in the magnetosphere with the 
initial twist ■00 — 0.2 and a discharge voltage V{u) — const. Time is expressed in units 
of ty — ji/cRV. Immediately after the starquake, the front emerges from the star at the 
magnetic equator and begins to expand outward and poleward. The front reaches m = (the 
magnetic axis) and erases all of the magnetospheric current in a finite time tend = ipoty. The 
dashed line shows the slope of 'U^(t) when -u^ approaches 0. This slope equals — (2-(/'o)^^- Upper 
panel: evolution of the dissipation power L (magnetospheric luminosity). L is expressed in 
units of Lv = V/ = cixV /AR^. The dashed curve corresponds to the dashed line in the bottom 
panel. Both panels show the model with ■^o = 0.2. Similar plots for twists with different 
initial ^00 = const are obtained by simple stretching of the time coordinate t — > t('0o/O.2); 
the evolution slows down by the factor ■0o/O.2 for stronger twists. 
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model, which assumes = 0. The drop appears in the numerical simulation at a finite 
M <^ 1 because the numerical model assumes a finite j^. The twist with ipolu) = const has 
the current density near the axis j{u) oc and hence, in a small polar region where j ^ j^,, 
$e must drop. In the limit jj, —>■ 0, this effect disappears in the sense that the polar region 
with j ~ jj, shrinks to one point u = 0. 

The twist behavior on the axis becomes important when the spindown of the star is of 
interest. The model in Figure 5 assumes a finite j^,, however, it does not take into account 
rotation of the star; therefore j(0) = and $e(0) = 0. Rotation with angular velocity 
Q implies the additional twisting and opening of the field lines that extend to the light 
cylinder R\c = c/Q. This persistent "external" twisting induces small but finite currents 
on the magnetic dipole axis, along the bundle of open field lines. The open bundle has the 
parameter uic = R/R\c ^ 1- In particular, magnetars have u\c ^ 10^"^. As long as the 
behavior of the closed magnetosphere u > u\c is concerned, the rotation can be neglected, 
and j,t — » is a good approximation. 

5.4. j-bundle with growing twist 

As the cavity expands, the current- carrying region u < U-^ shrinks. We will call the 
current- carrying bundle of magnetic field lines "j-bundle," for brevity. The numerical model 
of § 15.31 (Figs. 5 and 7) assumed that the discharge voltage is the same for all magnetospheric 
fiux surfaces, V{u) = const. Then the twist remains static inside the j-bundle. It freezes and 
waits while it is eaten by the expanding front 

In contrast, if V{u) ^ const, the twist in the j-bundle will change linearly with time as it 
waits for the front to come. This change is described by equation fl27j) [note that $e = 
inside the bundle u < u^] or equation fl39l) . The twist amplitude ip decreases if V < and 
grows if V > 0. Observational data (discussed below) suggest V > and the growth of ip 
near the axis u = 0. Despite the twist growth at small u, its total energy E^^ is decreasing 
with time as the j-bundle shrinks. This evolution is consistent with the energy conservation 
law (eq.EID. 

For illustration, we calculated the same model as in Figure 5 but with new threshold 
voltage V{u) = (0.04 + 2'u)^/^V, where V approximately equals the average of V{u) in the 
magnetosphere < m < 1. Figure 8 shows the evolution of the twist amplitude ip in this 
case and compares it with the case of V(u) = const. 

The twist evolution is described by simple analytical formulas when the j-bundle is 
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Fig. 8. — Evolution of the twist profile ip{u). (a) Model with V{u) = const (same model as 
in Fig. 5). (b) Similar model but with V{u) = (0.02 + 2w)^/^V. Solid curves show ijj{u, t) at 
three different moments of time t. Time is expressed in units of ty — /i/cRV in panel (a) 
and ty — /i/cRV in panel (b). Dashed curve shows ■^^ = i^iu*) — the twist amplitude at the 
boundary of the cavity throughout the entire history of its expansion from u^, — 1 to u^, — 0. 
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narrow, ^ 1. In the leading order of ^ 1 the front equation becomes 

-dr = - (4i^/c^)/;KKUl1„,i l^ + °<'">i- -*«^- 

where prime denotes the differentiation with respect to u. The relation between / and ijj 
(eq. [22]) gives Io{u) = {cfiipo/4:R'^)u'^ + 0{u^) where ipoi^) = ■ip{u,0) is the initial amplitude 
of the twist. Then we get, 

^ ^ "[2/iV;o(0)/cV+2V'(0)t + ^^'''^^^ • ^^^^ 
If V'(0) 7^ 0, integration of this equation yields 

u.it) = ^ In ^^^±^, to - (42) 
^ ' 2V'(0) to + t ci?V'(0)' ^ ' 

where tend is the time when the front reaches m = 0; it is finite in all cases. 

The twist amplitude 4>{u, t) inside the j-bundle grows according to equation (139|) (unless 
ip reaches ■i/'max ^1), 

^(0,t) =V^(0,0) (1 + ^ ) , 0<t<tend. (43) 



tr 



It grows by a large factor if tend ^ ^o- 



For example, consider a model with linear V{u) = V(0) + V'{0)u a.t u < u and constant 
V{u) = V{u) at u > u. Suppose V{u) ^ V(0). Then we find. 



V'(0,tend) 



2fM_i 
V(0) 



> 1. (44) 



The drop in V(mv,) at small delays the arrival of the front to m = and gives an exponen- 
tially longer time for the twist growth at u = 0. Thus, even a small twist with tpQ <^ 1 can 
grow to 'i/'max iuside the j-bundle. Further growth is impeded by the MHD instability. 



5.5. Localized starquakes 

Starquakes may rotate part of the crust and leave the rest of it untouched. Suppose 
that a ring U2 < u < ui has been rotated. Then the twist "0 7^ is created only in the 
region U2 < u < Ui. This implies I{u > ui) = 0, i.e. the net ejected current is zero, and 
hence the current density must change sign at some Um in the region U2 < u < Ui. The 
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Fig. 9. — Evolution of the magnetospheric twist that is created by rotation of a crustal 
ring. The initial twist amplitude ipQ ~ 0.2 in the region 0.3 < u < 0.7 and close to zero 
outside this region. The corresponding initial current distribution dlo/du is shown by dotted 
curve. Current is plotted in units of J = cfi/iR"^. The discharge voltage V{u) = (| + u)V 
is assumed in this model [V is the average of V{u) in the magnetosphere < m < 1]. Solid 
curve shows the current distribution dl/du at time t = O.OOSty {upper panel), t = 0.007tv 
{middle panel) and t = O.OSty {lower panel) , where ty = fi/cRV. 
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Fig. 10. — Upper panel: evolution of the dissipated power L in the model shown in Fig. 9. 
L is expressed in units of Ly = IV = cfiV/AR"^. Time is expressed in units of ty = fi/cRV. 
Lower panel: outward-propagating front that erases the ejected positive current (cf. Fig. 9). 
This front starts at Um ~ 0.63 where the ejected current density j oc dio/du changes sign. 
Vertical dotted hnes mark stages I-III in the twist evolution (see the text) . 
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current function Io{u) reaches a maximum at Um- An example of such a locahzed twist and 
its evolution are shown in Figure 9. Voltage V{u) = (| + u)V 7^ const is assumed in the 
model. 

Three stages may be noted in the evolution of the ring twist: 

I. One current front is launched from u = 1. It immediately jumps to the twist boundary 
Ml, and continues to propagate toward smaller u. In addition, two divergent current fronts 
are immediately launched from u = Um- Thus, two cavities form in the magnetosphere. Both 
cavities grow until they merge: the two fronts erasing the spike of negative current (Fig. 9) 
eventually meet and "annihilate." Stage I ends at this point (at time t ^ O.OOSty where 
ty = fi/cRV). 

II. The merged cavity continues to expand poleward and erase the remaining spike of 
positive current at smaller u. All of the initially injected current is erased (and stage II ends) 
at t ^ O.OStv 

III. The current front proceeds toward the axis, erasing the currents that have grown 
there (from zero) since the beginning of the twist evolution. 

Voltage V{u) 7^ const was chosen in the model to allow the twist growth near the axis. 
The simulation shows, however, that the growth is slow compared to the expansion of the 
cavity. Only at the last stage III does the grown current create a significant twist ijj near 
the axis and somewhat decelerate the expansion of the cavity. 

Figure 10 shows the evolution of the twist luminosity L{t). Its decrease is quickest 
during stage I, as the spike of negative current is quickly erased. Figure 10 also shows the 
propagation of the current front u^{t) that starts at Um and moves to smaller u, erasing the 
ejected positive current. 

Instead of a ring U2 < u < ui, an axisymmetric localized starquake may rotate a cap 
u < Ui- The cap-twist and ring- twist evolve in a similar way. In particular, stages I and II are 
similar. However, stage III is absent for the cap-twist. It has U2 = 0, i.e. the initially ejected 
positive currents occupy the enire region around the polar axis. Erasing these currents takes 
a longer time. As a result, the twist grows to a large amplitude at m <^ 1 before the cavity 
reaches the axis (§ 15. 4p . Then a maximally twisted narrow bundle with ip = V'max forms. 



- 30 - 



6. Observational effects 
6.1. Luminosity 

Energy dissipation in the untwisting magnetosphere is confined to the bundle of current- 
carrying field lines u < u^, (j-bundle). It can be a bright source of radiation, with luminosity 
equal to the rate of Ohmic dissipation L. The luminosity generally decreases as the magne- 
tosphere untwists. Examples of the evolution of L{t) are shown in Figures 7 and 10. 

A large fraction of the dissipated power may be radiated quasi-thermally at the foot- 
prints of the j-bundle as the accelerated magnetospheric particles run into the star (BT07), 
creating a hot spot on the surface 9 < 9-^. The area of this spot is given by, 

A ^ 7r{R sin 9^f = ttR^u^. (45) 

As the cavity expands in an untwisting magnetosphere (Fig. 6), the spot shrinks. 

The evolution of a narrow j-bundle (m^ <^ 1) with a uniform twist ip (e-g- ~ V'max^ 1) 
is described by simple formulas. The free energy of the twist (eq. [H]) is then given bjo 

Et. = £^l' du ^ ^ 4 X 10^^ Bl, Rl ^2 ul erg, (46) 

where Rq = R/10^ cm, = Bpoic/lO^'^ G, and -Bpoie = 2ji/R?. The luminosity of the 
j-bundle with ip ^ const can be immediately calculated for a given voltage V(u), 

The simplest model with V{u) = const gives 

L = Vh^ -^ipVul^l.3x 10^^ BuRQipVguleig s-\ (48) 
4ii^ 



where Vg = V/10^ V. The j-bundle with V{u) = const shrinks according to equation (138|) . 
which yields at Mj, ^ 1 

dUi, cRV , , cRV , , 

— ^ ~ ^ uAt) ^ tend - 1) . 49) 

dt ^ ' 2^i) V end ; \ > 



^"Only 1/4 of i?tw resides in the region u < 3/4 of the twist energy is contained in the potential region 
u> Ui,. The fact that V x B = in the potential region does not imply that = is determined by 



eq. ©. 
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This equation is also easy to derive from dEt^/dt = —L, using equations (H6|) and fHHj) . The 
evolution timescale of the luminosity is given by 

Approximate formulas can also be derived for V{u) ^ const. For example, for V{u) = 
V(0) + V{0)u we find 



1 + 



/2V'(0) 1 \ 2^ 
V3V(0) 2m J ^ 



(51) 



This equation again assumes ip^u < n^) ~ const. It may approximately describe e.g. the 
maximally twisted j-bundle, ip ~ V'max- 

The decay of the twist luminosity L{t) was previously estimated assuming that the 
current is decaying uniformly in the twisted region, which gave a linear L{t) oc t — tend 
(BT07). The electrodynamic theory developed in this paper shows that the untwisting is 
strongly non-uniform: the twist is erased by the propagating front that resembles a shock 
wave. The speed of this front depends on the initial twist configuration ipo{u). Similar to 
the simple estimate L{t) oc t — tend, we find that the twist is erased in a finite time and L{t) 
vanishes at tend (unless new starquakes occur). However, no universal linear shape of L(t) is 
predicted. In some cases L(t) may resemble a linear decay (e.g segment II in Fig. 10 is almost 
linear in a linear plot). In observed sources, L{t) was close to linear in AXP IE 1048.1-5937 
(Dib, Kaspi, & Gavriil 2008) and non-linear in XTE J1810-197(Gotthelf & Halpern 2007). 



6.2. Nonthermal radiation 

The energy released in the j-bundle can power nonthermal magnet ospheric emission. 
Such emission is observed in most magnetars. Two distinct non-thermal components are 
detected in their spectra: (1) soft X-ray tail that extends from 1 keV to ~ 10 — 20 keV with 
a photon index F ~ 2 — 4 (e.g. Woods & Thompson 2006 and refs. therein), and (2) hard 
X-ray component that extends to ~ 300 keV with F ~ 0.8 — 1.5 (e.g. Kuiper et al. 2008 and 
refs. therein). 

The 1-20 keV tail is usually explained by resonant scattering of thermal radiation by 
the magnetospheric plasma (Thompson et al. 2002; Lyutikov & Gavriil 2006; Fernandez & 
Thompson 2007; Nobili, Turolla, & Zane 2008; Rea et al. 2008). Ions resonantly scatter 
thermal photons near the star where heB/niiC ~ keV, and scatter at radii r ~ lOi? where 
heB/nieC ~ keV. The growth of cavity in the untwisting magnetosphere implies that the 
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scattering in the inner magnetosphere is suppressed, because the dense plasma is confined 
to the narrow j-bundle. At larger radii r ~ -R/u^, the j-bundle broadens and forms an outer 
corona that subtends a large solid angle as viewed from the star. The cyclotron energy in 
this region is 



and resonant scattering by can give 1-20 keV photons. The luminosity expected from a 
strongly twisted j-bundle with ~ 0.1 is consistent with the typical nonthermal luminosity 
of magnetars, L ~ 10'^^ erg/s (see eq. HH] and substitute the typical yU ~ 3 x 10'^^ G cm'^ 
and ip = 1 — 2). If no new starquakes occur, the j-bundle must shrink toward the magnetic 
axis with time, and will be reduced below 0.1. Then the resonant scattering must be 
suppressed. This suppression is caused by two reasons: heB/rrieC in the outer corona r ~ 
R/u^, decreases below keV (it is proportional to ul), and the power dissipated in the j-bundle 
becomes small as L oc m^. 

The j-bundle with u^, ~ 0.1 — 0.2 dissipates sufficient energy to explain also the hard 
X-ray component. This component was detected in three anomalous X-ray pulsars (AXPs) 
and two soft gamma-ray repeaters (SGRs) (see Kuiper et al. 2008 for a recent review). 
In AXPs, the hard X-ray emission has a huge pulsed fraction, approaching 100% at high 
energies (Kuiper et al. 2006; den Hartog et al. 2008a,b). This may be explained if the 
emission is produced in the narrow j-bundle near the star. 

Observations indicate that the nonthermal emission can be stable on timescales as long 
as a decade (den Hartog et al. 2008a). There may be two reasons for this stability: (1) 
The discharge voltage V is relatively low (below 1 GeV) and the j-bundle is relatively thick, 
~ 0.2. Then the untwisting timescale becomes long (see eq. [50]). (2) The j-bundle is 
kept in a quasi-steady, maximally twisted state ~ V^max by frequently repeating (possibly 
continual) shearing motion of the crust in a fixed region u < u^,. 

The plasma filling the j-bundle may also produce optical and infrared radiation by 
mechanisms discussed in BT07. Besides, it can be a bright source of radio waves, as suggested 
in section § 17.21 



The magnetospheric twist is expected to impact the spindown rate when the twist 
amplitude is large, ip ^ 1. This impact may occur in two ways. 

(1) The strong twist infiates the poloidal field lines and increases the magnetic field at 




(52) 



6.3. Outer magnetosphere and spindown rate of the star 
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the light cyhnder, which leads to stronger spindown torque acting on the star (Thompson 
et al. 2002). The poloidal inflation is common for twisted configurations. It is seen e.g. in 
the self-similarly twisted dipole (Wolfson 1995). A similar inflation must occur when the 
currents are confined to the j-bundle. Its calculation will require a full nonlinear model. Here 
we limit our consideration to simple estimates. 

One can think of poloidal inflation as an increase of magnetic dipole moment with 
radius. This effect is easiest to evaluate for moderate twists ip < 1. The current dl flowing 
along a bundle of twisted fleld lines du creates a toroidal current dl^j, ~ dl ip{u) /2t: . The 
fleld lines carrying this current extend to -Rmax = -R/m, and the dipole moment created by 
dl^ is du ~ dl^R^g^^/c. Integrating over m, we obtain the net change of dipole moment of 
the star due to the twist 

--^/^- 

The j-bundle u\c < u < with a uniform twist increases the dipole moment of the star 
by 

^~^ln^. (54) 

This effect is quadratic in and quickly becomes small for < 1. For strong twists ~ "^max, 
the estimate (IMI) must be replaced by a full nonlinear calculation. Qualitatively, it suggests 
that Afi/ fi is reduced as the j-bundle shrinks (u^ decreases). Therefore, the spindown torque 
is expected to be reduced with time and gradually come back to the standard dipole torque 
as Uic- 

(2) When the twist has grown to ipraax, it will begin to "boil over" through a repeated 
instability. The energy that would be stored in the magnetosphere if ip kept growing above 
V^max is then carried away by an intermittent magnetic outflow. The outflow may also carry 
away a signiflcant angular momentum. 

The outflow can be generated where the overtwisted field lines open up. This opening 
occurs because the twist is constantly pumped near the axis by V > 0. A possible structure 
of untwisting magnetospheres with V' > is schematically shown in Figure 11. It resembles 
the picture of a rotationally powered pulsar (see e.g. Arons 2008 for a review), however the 
opening is caused by the internal resistive dynamics of the magnetosphere and depends on 
the profile of V{u). By contrast, in ordinary pulsars the twist is pumped by the star rotation. 

Mechanisms (1) and (2) start immediately after the starquake if it implants a strong 
initial twist ipo > ^ near the magnetic axis, li ipo < 1, the spindown torque may not be 
affected until grows to ~ 1, which takes time (cf. eq. [27] or eq. [39!) 
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Fig. 11. — Possible structure of untwisting magnetospheres (poloidal cross section). The 
growing twist near the axis inflates the outer magnetosphere until its field lines open. Then 
the j-bundle {shaded in yellow) becomes confined between the last closed flux surface {blue 
curve) and the inner cavity {red curve), which expands with time. An equatorial outflow is 
expected to form just outside the last closed flux surface. The outflow is driven by dV/du > 0, 
which forces the twist to grow until part of the overtwisted fleld lines reconnect away from 
the star. The small circle in the center (shaded in cyan) shows the neutron star. 
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where V = dV/du is evaluated near the axis u ~ 0; V may be large, because u is small. For 
example, a change in V from V = 10^ V at u = to 2 x 10^ V at m = 10~^ corresponds to 
V' = 10^^ V. The delay is observed in some objects, with a characteristic tdeiay ~ 10^ s (e.g. 
Gavriil & Kaspi 2004). This requires V ~ 10^^ V. 

In contrast to luminosity L{t) (which always tends to decrease after the starquake), 
the torque behavior is generally non-monotonic. If ipo < the torque is expected to grow 
as ip grows in the shrinking j-bundle. (Such an anti-correlation between the torque and 
the X-ray luminosity was observed in IE 1048.1-5937, see Gavriil & Kaspi 2004). Once 
ijj has reached 

max ~ 1 the torque should start to decrease, because the shrinking of the 
j-bundle at constant ip ~ V'max leads to the reduction of A/x. The power of the outflow from 
the magnetosphere is also expected to decrease. The simultaneous decrease in torque and 
luminosity was observed, e.g., in XTE J1810-197 (Camilo et al. 2007). 

7. Untwisting magnetosphere in XTE J1810-197 

XTE J1810-197 is an anomalous X-ray pulsar with period P = 5.54 s and estimated 
dipole magnetic moment fi ~ 1.5 x 10^^ G cm^, which corresponds to the surface field at the 
polar cap 5 ~ 3 x 10^^ G (Gotthelf & Halpern 2007). An X-ray outburst was detected from 
this object in January 2003 (Ibrahim et al. 2004). Its luminosity approximately followed an 
exponential decay on a timescale of 233 days for 3 years (Gotthelf & Halpern 2007). During 
the X-ray decay, the source became radio-bright (Halpern et al. 2005) and switched on as 
a powerful radio pulsar with unusual spectrum and pulse-profile variations (Camilo et al. 
2007). The spindown rate of the star dramatically increased following the outburst. In the 
subsequent years, the object gradually evolved toward its quiescent (pre-outburst) state. 

These observations clearly indicate that the magnetosphere of XTE J1810-197 changed 
in the outburst, i.e. the footpoints of field lines must have moved, imparting a twist to the 
magnetosphere. The observational data give significant hints about the twist geometry and 
evolution: 

(1) The change in spindown rate suggests that the open field-line bundle was affected 
by a strong twist ip ^ 1 near the magnetic dipole axis. On the other hand, ip is limited to 
"^max = C^(l) by the MHD instability. Therefore, we infer ip = 0{1) near the axis. 

(2) Already one year after the outburst, the object luminosity was below 10^^ erg/s. 
The theoretically expected luminosity from a global twist with ip ^ 1 (eq. HHj) would be 
much higher: L ~ 3 x 10^^ Vg erg/s, and it would decay much slower than observed (eq. [50|l . 
Therefore, we conclude that the twisted region was small: the current-carrying field lines 
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formed a narrow bundle emerging from a small spot on the star surface. 

(3) Remarkably, the spot was discovered: a hot blackbody component with a small 
emission area was found in the X-ray spectrum following the outburst (Gotthelf & Halpern 
2007; Perna & Gotthelf 2008). Its emission area shrank with time until the spot became 
barely detectable (Fig. 12). 

(4) The X-ray and radio pulse profiles had almost simultaneous peaks, consistent with 
the X-ray hot spot being near the magnetic dipole axis (Camilo et al. 2007). 

These observations are consistent with the theory developed in the present paper. The 
hot spot is explained as the footprint of the j-bundle on the star (§ 16.11) . It has a sharp 
boundary and is shrinking with time as the j-bundle shrinks toward the magnetic dipole axis. 
The large amplitude of the twist ~ 1 may have been created initially by the starquake, but 
not necessarily: can naturally grow to ^ 1 following the starquake (§ 15.41) . A quantitative 
comparison of the model with the observational data is given below. 



7.1. X-ray emitting spot 

Useful preliminary estimates can be made if we assume a uniform voltage V{u) ~ const 
and uniform twist ip{u) ~ const across the j-bundle m < (§ 16. ip . Then the produced 
luminosity L{t) is given by equation PHI) . Suppose a large fraction of L is emitted thermally 
at the footprint of the j-bundle. The estimates for L (eq. HHl) and its evolution timescale tev 
(eq. ISUI) for a given spot area A (eq. HSl) may be compared with observations. For example, 
in the fall of 2004, the spot area was A ^ 10^^ cm (and hence ~ 0.03). One then finds 
that a strongly twisted j-bundle (t/; ~ 1 — 1.5) with voltage Vg ~ 3 explains both observed 
L «i 2 X 10^^ erg/s and tev ~ 0.6 yr0 

The model V(u) ^ const gives good estimates for L, A, and tev, however, it has a 
drawback: it is unable to describe the possible growth of from a smaller tpQ before the 
j-bundle became maximally twisted. The growth occurs if V > (§ 15. 4p . Therefore, we 
adopt a slightly more general model that includes the next (linear) term in the expansion of 
V(m) near u = 0, 

V(m)^Vo + V'm, m<1. (56) 
V is unknown and probably large near the axis (see the text after eq. ESD • 

A simplest model of a twisted magnetosphere has four parameters: magnetic dipole 



"We use here « 3 (Gotthelf & Halpern 2007). 
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Fig. 12. — Comparison of the model with the observed evolution of the area A and luminosity 
L of the hot spot formed after the outburst in XTE J1810-197. The data [open squares) 
are from Gotthelf & Halpern (2007). Dashed line shows the object luminosity in quiescence. 
Solid curve shows the theoretical model (see the text after eq. 1561) . 
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371) that we solve numerically. If ijj reaches 
The luminosity of a maximally twisted 



moment of the star /x, radius of the star R, the initial size of the j-bundle created by the 
starquake Uq = sin^ 6^0, and the initial amplitude of the twist ipo [it equals the angular 
displacement of the crustal cap rotated by the starquake; in our fiducial model, the starquake 
imparts a uniform twist ipoi'^ < ^o) = const]. The twist evolution after the starquake is 
controlled by two parameters Vo and V that specify voltage V(m) in equation (l56ll . The 
evolution is described by equations (l39l) and 
V'max = C'(l), it is assumed to stall at ipmay 
j-bundle, L{u^), is approximately given by equation (1^ . 

Figure 12 compares the model with observations of XTE J1810-197. The star is assumed 
to have = 1.5 x 10^^ G cm^ and R = 9 km. The parameters of the starquake are uq = 0.15 
and ipQ = 0.5. The discharge voltage is assumed to drop linearly from 5.5 GeV at u = 0.15 
to 1 GeV at u = (which corresponds to Vq = 10^ V and V = 3 x 10^° V). Although we 
did not attempt a formal fitting of the data, the figure suggests that the model is successful 
in explaining the evolution of L{t) and A{t). The data firmly constrain the voltage to be in 
1-6 GeV range, which is close to the theoretical estimate (eq. l33l) . 

The accuracy of the model is limited by several idealizing assumptions: uniform initial 
twist in the starquake region, the linear form of V{u), and axial symmetry. The actual value 
of fi is probably smaller than inferred from spindown measurements (§ 16.31) . The linear twist- 
evolution equation becomes approximate when ~ 1. Note also that our simplest fiducial 
model assumes that all energy dissipated in the j-bundle is emitted thermally at one (e.g. 
anode) footprint. A more realistic model can predict a more complicated spectrum with both 
footprints emitting. Finally, the hot-spot emission may be significantly anisotropic (Perna 
& Gotthelf 2008), which may increase or reduce its apparent luminosity depending of the 
average inclination of the spot to the line of sight. 

Note that the observed spectrum was well fitted by a two-temperature blackbody. The 
second blackbody component was cooler and had emission area comparable to that of the 
star. No non-thermal component from magnetospheric scattering was required by the data. 
The weakness of the scattered component may be explained by the small size of the 
j-bundle: the star's radiation does not scatter in the cavity around the j-bundle, and the 
probability of scattering in the narrow bundle is low 



^^In the simulations, we set 'i/'max = 1-5. The exact value of the stability threshold needs to be calculated 
in the full nonlinear twist model. 



- 39 - 



7.2. Radio pulsations 



In ordinary radio pulsars, radio emission is believed to be produced by the open field 
lines passing through the light cylinder, because they are the only part of the magnetosphere 
that carries electric currents. In magnetars, these field lines have a tiny uic = R/R\c ~ 10~^ 
where R\c = c/Qis the light- cylinder radius for a star rotating with angular velocity Q. Radio 
emission from the open bundles of magnetars may be undetectable for two reasons: (1) 
discharge has a low threshold in magnetars, V ~ 10^ V. When this voltage is multiplied 
with the current in the bundle, lie ~ lufc = cfi/iRy^, one obtains the dissipated power 
Lie = V/ic ~ 10^^ Vg erg/s. This power appears to be too small to feed the observed radio 
luminosity of XTE J1810-197, Lj-adio ~ 10^° erg/s. (2) The radio beam may be narrow 
(because of small uic). Then the probability of its passing through our line of sight is small. 
This suggests that radio pulsations are hardly detectable when the magnetar is in quiescence, 
i.e. when its magnetosphere is untwisted. 

In contrast, after the starquake, the j-bundle forms. It is much thicker and more ener- 
getic than the bundle passing through the light cylinder and can produce much brighter radio 
emission with a much broader pulse. The untwisting magnetosphere in XTE J1810-197 had 
U-k/uic ^ 3 X 10^. The net current flowing in the j-bundle is ~ 10^ times larger than lie, 



This gives L/Li^ ~ 10^ (assuming a comparable discharge voltage at and uic). A small 
fraction eradio of this luminosity may escape as radio waves; eradio ~ 10~^ is consistent with 
observations. 

Radio waves are efficiently absorbed by the magnetospheric plasma. Only waves pro- 
duced close to the magnetic axis, u < Uesc, are likely to escape, while waves emitted at 
u > Uesc are trapped in the magnetosphere. The value of Uesc depends on the frequency of 
the wave and the state of the plasma (particle density and velocity distribution); Uesc might 
be inferred from the observed opening angle of the radio beam. The radio luminosity should 
quickly decrease when the j-bundle shrinks to < Wesc- 

Radio pulsar XTE J1810-197 is distinguished from ordinary radio pulsars by its hard 
spectrum, very strong linear polarization, and variable pulse profile. If its emission is pro- 
duced on the bundle of closed field lines with u ^ uic, it may be expected to be different from 
ordinary radio pulsars. The spectrum of radio waves may form as the sum of emissions from 
different (frequency-dependent) Uesc inside the j-bundle. The sporadic changes in the pulse 
profile may be caused by the instabilities in the maximally-twisted outer magnetosphere 




(57) 



(§[631). 
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Radio observations provided accurate measurements of the spindown torque acting on 
the star (Camilo et al. 2007). The torque was decreasing together with the X-ray luminosity 
2-3 years after the starquake. Its history at earher times was not observed; the torque is 
beheved to have increased after the starquake (Camilo et al. 2007). Such a non-monotonic 
evolution of the torque would be consistent with the theoretical expectations (§ 16. 3p . 

8. Discussion 

Electrodynamics of untwisting may be summarized as follows: 

1. — The twist evolution following a starquake is not diffusive spreading that was 
pictured previously. Instead, the twist current is sucked into the star, a cavity immediately 
forms in the inner magnetosphere and grows until it erases all of the twist (Fig. 6). A sharp, 
step-like drop in current density j is maintained at the boundary of the cavity. It is caused 
by the threshold nature of the discharge that conducts magnetospheric currents. 

2. — As the cavity expands and the j-bundle shrinks toward the magnetic axis, the 
twist amplitude ip in the j-bundle can grow. The growth occurs if dV/du > 0, i.e. the 
discharge voltage is smaller on field lines extending farther from the star. This is plausible 
if the current is conducted through dischargeo The j-bundle with the growing twist is 
shrinking faster with time, so that the total twist energy Et^^ decreases consistently with the 
Ohmic dissipation rate. 

3. — The growing twist that has reached the threshold for instability "i/'max = C^(l) 
is expected to "boil over" and drive an intermittent outflow of magnetic energy from the 
star. The twist in the j-bundle then remains near iprn&x for the rest of its lifetime. It may 
be regulated by the limit-cycle instability — the repeated growth of ip to ipma^ followed by 
a sudden reduction of ip below ipmax- The value of ipma.x and the nonlinear behavior of the 
outer magnetosphere at ~ V'max needs to be studied further using numerical simulations. 

4. — In addition to rare large starquakes, the magnetosphere can be gradually twisted 
by the continued motion of its footpoints, either plastic or through a sequence of small 
starquakes. The continued footpoint motion leads to either a very strong twist or a tiny 
(negligible) twist that has no observational effects. A quasi-steady state with a non-negligible 
twist amplitude is possible only with ip ~ V'max- 



^•^In contrast, if the current-carrying charges were lifted from the star by voltage ([32|l . the voltage would 
be larger for field lines that extend to larger altitudes, i.e. dV/du < 0. Then the twist on the j-bundlc would 
diminish with time rather than grow. 
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Recent observations of XTE J1810-197 are particularly useful for testing the untwisting 
theory, for a few reasons: (1) XTE J1810-197 displayed a clean post-starquake evolution, 
which was observed for years uninterrupted by new starquakes. (2) The low level of quiescent 
luminosity from this object allows one to see clearly the shrinking hot spot on the star, and 
its evolving luminosity was tracked from 10^^ erg/s down to 10^^ erg/s. (3) The detection of 
radio pulsations and detailed measurements of spindown torque make this object yet more 
interesting for testing theoretical models. For these reasons, this paper focused mainly on 
XTE J1810-197(§ 7). Recently, an outburst was detected in the similar radio magnetar 
IE 1547.0-5408 (Camilo et al. 2008; Halpern et al. 2008). Its behavior appears to be more 
complicated than that of XTE J1810-197; apparently, episodes of repeated (and overlapping 
in time) activity occurred a few months after the outburst. The analysis of this object is 
deferred to a future work. 

Other, more active, AXPs and SGRs display a diverse and complicated behavior of 
X-ray luminosity, pulse profile, and spindown rate (see Woods & Thompson 2006; Kaspi 
2007; Mereghetti 2008 for reviews). Repeated starquakes of various amplitudes and possible 
plastic deformations of the crust make these objects more difficult to analyze. The continuing 
injection of a magnetospheric twist can slow down the decay of its luminosity. It can also 
affect the star's spindown in a more complicated way than described in § 16.31 for an isolated 
starquake. In general, twist injection should lead to higher X-ray activity and faster spindown 
(Thompson et al. 2002). Such a general correlation exists in the magnetar population 
(Marsden & White 2001) but not always observed in individual objects. Note that spindown 
is controlled by the X-ray-dim bundle of open field lines. The impact of a starquake on 
this narrow bundle may be immediate, occur with a delay, or never occur, depending on the 
starquake geometry and amplitude (cf. the case of a "ring" starquake in § 15.51) . Repeated 
starquakes may lead to a non-trivial relation between spindown and X-ray emission. 

Despite the diverse behavior of active magnetars, some general features may be inferred. 
It is clear that the observed sources do not have strong global twists, for two reasons. First, 
the evolution timescale of such twists would be too long, t 

ev 

(10 - lO'^)Vg^ yr, where 

Vg is the discharge voltage in units of 10^ V (see eq. |50]). In contrast, the magnetospheres 
of observed magnetars usually evolve on timescales ~ 1 yr or even shorter!^ Second, the 
luminosity produced by a strong global twist would be too high, L ~ 10^^ Vg erg s~^ (see 
eq. 148!) . It is 2 orders of magnitude higher than the typical observed luminosities of magne- 



^"'The theoretical iov would be reduced for a larger discharge voltage V ^ 10^ V, however it appears 
impossible to sustain such a voltage as it leads to runaway creation (BT07). Even if V ^ 10^ V 
were theoretically possible, strong global twists would still be ruled out by observations because they are 
over luminous. 
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tars, L ~ 10^^ erg/s (e.g. Durant & van Kerkwijk 2006). This leaves two possibilities: the 
twist is weak [ip oi localized to a narrow bundle of field lines. The changing spindown 
rate suggests a strong twist, at least near the magnetic dipole axis. Thus, observations ap- 
pear to support the picture of a strongly twisted j-bundle near the dipole axis. It is certainly 
supported by the observations of XTE J1810-197 (§ [7]). 

The localized strong twist may be created by starquakes localized to the polar region. It 
also tends to form dynamically from a weak global twist as the j-bundle shrinks to the axis. 
The luminosity and evolution timescale of an untwisting magnetosphere with u^, = sin^ ~ 
0.1 and ?/' ~ 1 is consistent with typical L ~ 10^^ erg/s and tev ~ 1 yr of active magnetars 
(see eqs. HHl and [50]) . Both nonthermal X-ray components in the magnetar spectra, 1-20 keV 
and 20-300 keV, can be produced by the j-bundle. They are likely emitted at different radii. 
If the hard component is produced near the star, where the j-bundle is narrow, a relatively 
narrow 20-300 keV pulse is expected. 

The resonant-scattering model for 1-20 keV radiation is consistent with the picture of 
a narrow j-bundle near the dipole axis. The cyclotron resonance for with keV photons 
takes place at radii r ~ 10-R. This means that the resonant scattering is confined to the 
field-line bundle with u = R/R 

max ~ 0.1 (§ 16. 2p . The luminosity of upscattered radiation 
Lsc is supplied by Ohmic dissipation of electric currents in this bundle. 

The footprints of a bundle with u ~ 0.1 form ~ 3 km spot on the star. A significant 
fraction of energy released in the j-bundle may be transported to its footprints and emitted 
there quasi-thermally. If the spot radiates a significant part of the bundle luminosity, L ~ 
-^q35 erg/s, then it must have a temperature kT ^ 1 keV. Such hot spots were observed in 
XTE J1810-197 and IE 1547.0-5408, but were not reported for other, brighter magnetars. 
The typical reported temperatures of blackbody components are 0.3 — 0.6 keV (e.g. Perna 
et al. 2001). The disappearance of hot spots in the presence of large Lsc may be caused by 
the outward drag created by resonant scattering, which suppresses the transport of released 
energy to the footprints of the j-bundle. 
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Appendix A: Twisted dipole magnetosphere 

Consider a weakly twisted dipole field B = Bq + B^. Its poloidal component Bq is 
that of untwisted dipole, 

5. = Be = (58) 
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The poloidal flux function is given by 

r'^^ sin^ 9 

f{r,d)= / / Brh^hed9d4> = 2nfi , (59) 

Jo Jo 

wliere = r sin 6 and hg = r. Since / = const on any flux surface, 

r 

= const = i?rnax, (60) 

sm 6 

which is the radius where the flux surface crosses the equatorial plane (the maximum radius 
reached by the flux surface). 

The magnitude of the twisted magnetic field is given by (neglecting terms C[i?^/i?^]), 

B{r, 6) ^ Bo = + 3 cos^ 9. (61) 

The twist angle if) (eq. [3]) is easy to calculate using 9 as the parameter along the field line 
and substituting dl = {B/Bg)d9, 

■,^r^'-f. (62) 

Jeo 

Here is the polar angle of the northern footpoint of the field line, and vr — is the polar 
angle of the southern footpoint. Substituting B^ = 21 /ch^ (eq. H]), Bq from equation fl58|) . 
and using equation fIBU]) . one finds 

4/ 

^ = — ^maxCOS^O, (63) 



where 



cos2eo=l-^ = l-y. (64) 

When the poloidal current / and twist angle if) are viewed as functions of m = and time 
t, the relation between ip and / (eq. [63]) becomes equation fl26|) . The free energy of a twist 
with given I{u) is found using equations f|T9|) and fl26|) . 



^tw(t) = — / l\u,t)Jl -. (65) 

C ./q V Ur U 



We also calculate here the integral that is needed in § 5.1, 



Bdl^ Bodl. (66) 

r>R Jr>R 
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It is taken along the field line outside the star. Consider the contour that closes the path of 
integration between the two footpoints along the surface of the star (and across the magnetic 
field) . The line integral of Bg along this closed contour vanishes (as follows from the Stokes' 
theorem and V x Bq = 0). This gives 

/ Bodl= BeRde = ^^i^. (67) 

Jr>R Jei ^ 

Here 6i and vr — 6'i are the polar angles of the field-line footprints on the star surface, and 

cos 6*1 = yjl — U. 



Appendix B: Derivation of the front equation 



The front is located on the flux surface u^, where ^eiu) deviates from V(m), current j 
jumps to ~ (Fig. 5). We consider below the limit — >• 0. Then the jump of j is described 
by the Heaviside step function 0(mj, — u) and j = in the region u^, < u < 1. 

This implies that I{u) = const for U-^ < u < 1, and this fact can be used to derive 
the dynamical equation for m*. Let us differentiate the twist evolution equation (l28l) with 
respect to u. Then the left-hand side vanishes in the region u-^ < u < 1 and we get 



d 



du \aR^I - u/uc du 



0, 



which implies 



(68) 



(69) 



du ^ 

where K is constant in the region u^, < u < 1. Integrating this equation and using the 
boundary condition $e(l) = we find the shape of $e(^); 



V{u) 
K 



X 1,1 

- + - In - 

X 2 1 



X 



Q < U < Ui, 
Ui, < u < 1 



(70) 



' J X=u/Uc 

The explicit expression for K can be found from the condition $e(M*) = V(u^ 



Uc 




UZ 



a/1 - U^/Uc 



'I 



,-1 



1 + Jl-U^/Uc 



(71) 
(72) 
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The numerical factor ^(u*) is well approximated by a simpler formula, 

1 + lAu^ 



1 



(73) 



The accuracy of this approximation is better than 2.2% for < < 0.65. The approx- 
imation is worst (30% accuracy) when = 1 and quickly becomes excellent as the front 
propagates away from the boundary. 

The twist evolution equation fl28l) can now be written as 



m 

dt 



4:Rjl - U/Uc 

4R 



< M < Mi 



Ui, < U < 1 



(74) 



where V = dV/du. Note that 81/ dt does not vary with time at m < m^, and hence I{u,t) 
can be obtained by simple integration. Thus we find the current function of the twist. 



Iiu,t) 



m 



cW{u) 



:t < U < U^, 



< U < I 



(75) 



where Io{u) = l{u,0) is the initial current function, and /^(t) = I[Ui,it),t] is the net current 
that flows through the magnetosphere at time t. Using the relation between / and ijj (eq. [26]) . 
we obtain equation (|39|1 . It gives a complete solution to the problem of twist evolution if we 
find Ui,(t). 

The equation for u^{t) can be obtained if one considers the evolution of I^. Its time 
derivative is given by 

dl^, dl . , dl . , dUi, 
which can be evaluated separately on each side of the front. 
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dh 
dt 



U^,-0 



dh 
dt 



dio 
du 



w*+0 



+ 



d cW(u 



du AR^/l - u/uc 



dUi, 



c^K 
4R ■ 



(77) 



(7J 



The two results must match because they describe the same /*(t). The jump in dl /du 
across the front (which corresponds to the jump in j) is compensated by the jump in dl /dt. 
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Equating the two expressions for c?/*/ dt, we find 



~dt 



a/1 - M^/Mc 



AR dio 
du 



+ 



d u^V'iu) 



du y/l - u/ 



(79) 



This ordinary differential equation describes the propagation of the front. Substitution of 
equation ( ITTi) for K{u^,) gives equation ( 1371) . 
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